Evaporative depolarization and spin transport in a unitary trapped Fermi gas 



o 
o 

(N 
o 

Q 



I 

a ■ 

■ 



I 

o 



(N 
> 

m 
vn 

rn 
O 

ON 

o 



■ 

X 



Meera M. Parish^^^ and David A. Huse^ 

'Department of Physics, Princeton University, Princeton, NJ 08544 
^Princeton Center for Theoretical Science, Princeton University, Princeton, NJ 08544 

(Dated: December 17, 2009) 

We consider a partially spin-polarized atomic Fermi gas in a high-aspect-ratio trap, with a flux of 
predominantly spin-up atoms exiting the center of the trap. We argue that such a scenario can be 
produced by evaporative cooling, and we find that it can result in a substantially non-equilibrium 
polarization pattern for typical experimental parameters. We offer this as a possible explanation for 
the quantitative discrepancies in recent experiments on spin-imbalanced unitary Fermi gases. 



I. INTRODUCTION 

Two-component atomic Fermi gases provide an ideal 
experimental system in which to investigate fermion pair- 
ing and superfluidity in a controllable manner.— For ex- 
ample, one can use a magnetically-tunable Feshbach res- 
onance to access the unitary regime, where the scatter- 
ing length diverges and one has a strongly-interacting 
fermionic superfluid that is 'universal'i^ Of particular in- 
terest is the case where there is a spin imbalance that 
frustrates pairing between fermion species, because this 
is a situation that arises in many fields of physics, ranging 
from QCD to superconductivityi^ Here, a central ques- 
tion has been: what is the critical spin polarization Sc at 
which pairing and superfluidity are destroyed for a uni- 
tary trapped Fermi gas at equilibrium? However, cur- 
rent experiments produce different answers. The exper- 
imental group at MIT finds that Sc ^ 77%^ii^^ while 
experiments on highly-elongated trapped gases at Rice 
University^ suggest that 6c is at least 90%. Moreover, 
even though the critical polarization is known to be a 
strong function of temperature, where Sc decreases with 
increasing temperature, the different (5c 's observed in ex- 
periment are unlikely to be caused by differences in tem- 
perature since both experimental groups have claimed 
that their temperatures are low enough to yield ScS that 
are close to the zero-temperature result. Thus, there is 
a real discrepancy in the measured Sc and a resolution 
of this problem has potentially important implications 
for the nature of the paired superfluid phase in a finite- 
sized system. Here we propose that the high Sc observed 
in the Rice experiment was due to their trapped spin- 
imbalanced gas being out of equilibrium. We show that 
a combination of the trap geometry and the evaporative 
cooling scheme implemented in the Rice experiment can 
induce a spin current along the trap axis, which in turn 
creates a substantially non-equilibrium polarization pat- 
tern that favors a superfluid at the trap center. 

To understand how such a spin current can be gener- 
ated, one must first examine the evaporative cooling pro- 
cess. Here, the temperature, and entropy per atom, of 
a trapped gas is lowered when the most energetic atoms 
escape over the "lip" of the trap — the route of escape 
with the lowest potential barrier to be surmounted. For a 
partially-polarized Fermi gas at temperature T, the rate 
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FIG. 1: (Color online) Schematic diagram of atom trans- 
port during the evaporative cooling process in trapped, spin- 
imbalanced, highly-degenerate Fermi gases in the quasi-lD (a) 
and 3D (b) regimes. The flux of spin-up atoms jo is exiting 
the gas over the lip of the trap that is located at the axial cen- 
ter {z = 0). This flux must be drawn from the fully-polarized 
(FP) normal region. Thus, for a quasi- ID gas (a), a spin 
current must flow through the intervening partially-polarized 
(PP) normal region and across the normal/superfluid(SF) in- 
terface. In the 3D regime (b), the atoms can evaporate di- 
rectly from the surrounding fully-polarized layer. 



of thermal activation over this barrier is larger for the 
majority species by a factor of exp [(/X| — ^i)/kBT], as- 
suming the two species are subject to the same trapping 
potential, where /i| are the chemical potentials for 
the majority and minority species, respectively. Thus, at 
low T, the flux of evaporating atoms passing over the lip 
is essentially fully polarized, and we have evaporative de- 
polarization in addition to evaporative cooling (as stated 
in Ref. jj). The Rice experiments we are considering^ 
had both a long, thin, high-aspect-ratio optical trap and 
a nonuniform magnetic field that contributed to the axial 
confinement of the gas. As a result, the lowest barrier for 
the atoms to escape from the trap was at the axial cen- 
ter, with the atoms escaping over this "lip" in the radial 
direction (and downwards, due to gravity; see Fig. [Ija)). 
Elsewhere in the trap the barrier to escape was signifi- 
cantly higher, so at their lowest temperatures, essentially 
all of the evaporating atoms escaped radially at the axial 
center of their high-aspect-ratio trap. 

To achieve low temperatures in the Rice experiment, 
the height of this barrier was lowered by reducing the in- 
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tensity of the optical trap. It was at the lowest tempera- 
tures that the unusually large Sc was observed, along with 
strong deviations from the equilibrium local density ap- 
proximation (LDA) in the shapes of the regions occupied 
by the superfluid and partially-polarized normal phases.® 
What we propose happened here is that the evapora- 
tion, with the t atoms rapidly escaping radially over the 
trap lip, greatly depleted any excess unpaired f atoms 
from the axially central region of the cloud occupied by 
the paired superfluid phase. This depletion, which is ap- 
parent in the in situ density measurements;^'^ substan- 
tially suppressed (//| ~ Mi) that region (evaporative 
depolarization). The flux of evaporating "f atoms over 
the lip then had to come from the fully-polarized nor- 
mal regions at the axial ends of the cloud and be driven 
through the partially-polarized region and across the nor- 
mal/superfluid interface by a substantial axial gradient 
of (/^i — fii). This resulted in the partially-polarized re- 
gion of the cloud being much smaller in axial extent than 
it would be at equilibrium, which is the strong deviation 
from equilibrium LDA that we will focus on in this paper. 
There was another important deviation from equilibrium 
LDA: the aspect ratio of the central superfluid region was 
substantially reduced from that of the cloud as a whole; 
it is this latter feature that was emphasized in Ref. H. 

We emphasize that the non-equilibrium scenario de- 
scribed above assumes that the atom transport is effec- 
tively one-dimensional (ID) along the axial (z) direction, 
as illustrated in Fig. [IJa).^° Thus, it is only appropriate 
for high-aspect-ratio gas clouds with sufhciently low par- 
ticle number, like those in the Rice experiments (aspect 
ratio > 30 : 1) at the lowest temperatures. For higher 
temperatures (less evaporative cooling), where the parti- 
cle numbers are larger, there is a fully-polarized layer of 
t atoms fully surrounding the cloud (Fig. [Tfb)), as can 
be seen in in situ density measurements In this case, 
the evaporation will simply draw | atoms from this fully- 
polarized layer rather than driving a spin current through 
the superfluid and partially-polarized normal regions, so 
a strong chemical potential gradient is not produced. The 
MIT experiments, on the other hand, used trapped gases 
with a much lower aspect ratio (^5:1) and at least 
an order of magnitude more atoms in the cloud. Thus, 
the 3D scenario for highly-polarized gases depicted in 
Fig. [Ijb) remains correct for these experiments, even at 
the lowest temperatures, and so their observed 6c agrees 
with the equilibrium result from quantum Monte Carlo 
(QMC) calculations^^ 



II. MODEL OF EVAPORATION 

To model the spin transport in the Rice experiment 
at the lowest temperatures, we will assume that atoms 
are only removed from the gas at z = 0, the position 
of the trap lip, and that only f atoms are evaporating 
(Fig. ma)). Also, we approximate any nonuniformities 
in the local chemical potentials as ID, so that we only 



consider non-equilibrium axial gradients in fi^. Since 
we are primarily interested in the non-equilibrium polar- 
ization pattern within a trap, and this is not strongly 
dependent on the local temperature, we will ignore any 
gradients in temperature. We further assume that the 
atomic gas cloud is in mechanical equilibrium through- 
out the evaporation process, since the time required to 
equalize a pressure difference in the cloud is much smaller 
than both the duration of the evaporation and the time 
for spin diffusion along the length of the cloud. This 
implies that the local pressure P always satisfies: 
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where n = n-^ + ni is the local total density and V{z) is 
the axial trapping potential. 



III. SPIN TRANSPORT IN A UNITARY 
NORMAL GAS 

We can determine 5c of a spin-imbalanced Fermi gas in 
the presence of a spin current by focussing on the unitary 
normal gas at low temperatures. Here, the equation of 
state for the pressure at T = is accurately known from 
QMC calculationsiiii^ and is given by 
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where m is the atomic mass, m* ~ 1.09to ~ to, A ~ 0.99 
and F ~ 0.14. For the spin transport that we are in- 
terested in, we also expect Eq. ^ to provide a rea- 
sonable approximation for the pressure at low temper- 
atures within the degenerate regime T < sp] /ks = Tp-\, 
where £F,a = fi.^(67r^ncr)^/^/2TO. In addition to Eq. ([T]), 
we require equations for the spin density — ~ 
and spin current density jg = j-^ — j| . We will assume 
that the "DC" spin transport in the partially-polarized 
normal Fermi gas is diffusive, with the dissipation being 
driven by the interspecies interactions. Indeed, this situ- 
ation may be regarded as the cold-atom analogue of spin 
Coulomb drag in electron systems By transforming to 
the inertial reference frame where locally j' = = 0, 

i.e. where there is no net transport of the total mass den- 
sity and the motion is purely diffusive, we obtain current 
densities: 



"I dz dz 



(3) 



where —gf- is the local equilibrium density gradient and 
Dc is the diffusion constant for each spin. One can de- 
termine as a function of the local densities using the 
LDA equilibrium condition for the chemical potentials of 
each spin: 

djia_ _ dV 
dz dz 



(4) 
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Note that this automatically implies that there are no 
gradients in the chemical potential difference — /i^) 
at equilibrium. Combining Eq. ([3]) with the expression 
for the spin current density in this frame 



j's —j\ j'l — 



dz 



yields an expression for the spin diffusion constant: 
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Transforming back to the lab reference frame then gives 
us the spin current density: 



dz 



dz 



(7) 



where the average net velocity v = j /n. 

For a dilute gas with s-wave scattering, we have diffu- 
sion constants D^j = ^v'^t^, where v^- is the average ve- 
locity of the random particle motion and I/to- is the scat- 
tering rate of each species. In the regime T < Tpi < Tp^, 
the scattering rate for the J, atoms at unitarity isi^ 



1 ^Att^A^ m*{kBT)^ 



(8) 
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FIG. 2: (Color online) Total critical polarization 5c of a 
trapped quasi- ID spin-imbalanced Fermi gas as a function 
of the dimensionless spin current Q. The figure insets de- 
pict cross sections of the trapped normal gas at 5c for (5 = 
and Q = 1. In the latter case, the partially-polarized (PP) 
normal core has been distorted and shrunk compared to the 
fully-polarized (FP) region, in clear violation of LDA. Note 
that the depicted trap aspect ratio is at least 5 times smaller 
than that in the Rice experiment, and thus the apparent FP 
layer surrounding the PP core near the axial center {z — 0) is 
extremely thin, being only a couple of atoms thick, and does 
not violate our approximation of ID transport. 



By knowing T| and Va- , we can also extract an expression 
for the t scattering rate since the mean free paths = 
VrjTa are simply related by the densities: Z| = Z|(nf/n|). 
One can see this from the fact that = l/(Tcs?^-<T in 
a dilute gas, where the scattering cross section CTcs is a 
universal function of the density at unitarity. Also, we 
can approximate the velocities as Va ~ KkFa/fn- Thus, 
we obtain the 'universal' spin diffusion constant: 
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Note that the dependence of Dg on temperature is that 
of a Fermi liquid, as expected. Experiments to explic- 
itly measure this 'universal' transport constant would be 
welcome. 

Referring to Fig. [TJa) , if the fully-polarized normal re- 
gions at the axial ends are sufficiently large (i.e. if the gas 
has a sufficiently high global polarization), then we can 
treat them as stationary spin current sources and approx- 
imate the spin current through the partially-polarized re- 
gion as steady-state: 



dt 



dja_ 

dz 
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(10) 



In this case, by solving Eq. ([7]) with jg = j = Jq, together 
with Eq. ^ , we find that the deviation of the total polar- 
ization from equilibrium for a given local polarization at 



the trap center is determined by the simple dimensionless 
quantity: 
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where L is the total axial length of the cloud. 

The critical total polarization Sc for a trapped cloud 
is obtained by setting the local densities nj^/n-f = 0.44 
at the trap center, which is its value at the superfluid- 
normal transition in the uniform systemi^i This corre- 
sponds to the situation where a trapped normal gas is 
on the verge of forming a superfluid core. To obtain a 
quantitative estimate of Sc as a function of Q in a high- 
aspect-ratio harmonic trap, we determine the chemical 
potentials as a function of z and then we use LDA in 
the radial direction r to now include the radial trapping 
potential V±{r) and obtain the densities of each spin: 

ng{z,r)=ng[^i^{z)-V^{r),^i^{z)-V^{r)] . (12) 

Integrating these densities numerically then yields the 
total polarization S = {N-^ - Ni)/{Ni + N^). 

As depicted in Fig. [21 the critical polarization 6c cor- 
responds (by construction) to the equilibrium QMC re- 
suli^ii when (5 = 0, but it dramatically increases with 
increasing Q and asymptotically approaches 100% po- 
larization as Q —> oo. This increase in Sc is due to the 
partially-polarized region being compressed along the ax- 
ial direction relative to its extent in equilibrium LDA 
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(inset of Fig. [2]), an effect which has been observed in 
the Rice experiments. To obtain 5c — 90%, we re- 
quire an evaporation rate of order 10® atoms/sec if we 
use typical parameter values in the Rice experiment to- 
wards the end of evaporation: T jTp^ ~ 0.05, L ~ 1mm, 

rt| ~ lO^^cm^'^, h/m lO^/im'^/s, and a trap radius 
of 10 ^m. This appears consistent with the data in 
Ref. iSi since the actual evaporation rate is estimated to 
be roughly 10® atoms/sec at the time when the optical 
absorption images are takenJ^ 
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IV. SUPERFLUID/NORMAL INTERFACE 

For completeness, we now examine the effect of a su- 
perfluid core on the spin current in the trapped gas when 
5 < 5c- For this purpose, we consider spin transport 
through the superfluid/normal interface. We approxi- 
mate the interface to be a ID sharp step-function in the 
local superfluid gap A(z) and densities (Fig. [3]), which is 
reasonable when T is well below the temperature at the 
tricritical point Clearly, our assumption of mechani- 
cal equilibrium implies that the pressures in each phase 
are equal at the interface: Pm = Psf- Since we are fo- 
cussing on low r, we will use the T = expression for 
the superfluid pressure: 



FIG. 3: (Color online) Interface between the unpolarized 
superfluid and partially-polarized normal state at unitarity, 
where the local chemical potentials are shown for the inter- 
face at local equilibrium, with superfluid chemical potential 
2/isF = + Mj,^^'- The shaded energy range corresponds 

to the quasiparticle gap (the spin gap) of the superfluid, with 
A ~ 1.2^sF taken from T = QMC calculations.— 



current density of ] atoms flowing into the superfluid as: 



3] 



(27r)2fi 



dE / dk 



(14) 



2 (2mV'^ 5/2 

where ^ ~ 0.42 according to QMC calculations,^^ and the 
superfluid chemical potential 2^,sf = + Mj^^^''- We 

further assume that any quasiparticles in the superfluid 
region are rapidly evaporated so that we can neglect any 
quasiparticles incident on the interface from the super- 
fluid side. Effectively, this amounts to assuming a strong 
drop in the local temperature of the quasiparticles as one 
crosses the interface from the normal to the superfluid 
side. While this is not essential to our picture, this gives 
a simplification and likely captures what happens when 
the lip of the trap is held very low. 

The scattering problem at the superfluid/normal in- 
terface has already been examined in Ref. [l9i . However, 
they only focussed on thermal transport via quasiparticle 
transmission through the interface. In our case, we must 
consider both mass and spin transport across the inter- 
face, since they can both change the local polarization. 
We must therefore take account of any mass transported 
via Andreev reflection. For example, if the chemical po- 
tentials ^'^^\ M^^' ill the normal phase lie within the 
gapped region of the superfluid, like in Fig. [31 the An- 
dreev process results in a flux of pairs (mass) into the 

superfluid when 2^sf < + /^[^^ ^'^'^ ^ "-"^ 
the superfluid when 2^sf > + even at T = 0. 
Following the approach of Ref. [2^ we can write the 



where fc| (fc|) is the momentum normal to the interface 
of the incident f atom (reflected [ hole), E is the en- 
ergy and f{x) is the Fermi function. The coefficients 
v4| and correspond to the probability of Andreev re- 
flection and ordinary reflection, respectively. Clearly, 
when \E — ^sf\ < A, there is no quasiparticle trans- 
mission and we must have k-^ = k^B^ + ^i^f Also, 
for E < A + fj,sF (appropriate for low T), the current 
density of | holes crossing the interface into the super- 
fluid has the same expression as in Eq. p^ . but with 
1 and I interchanged, and /(x) f{—x)- Note that 
although our calculation of reflection probabilities relies 
on the mean-field Bogoliubov-de-Gennes equations, the 
chemical potentials, pressure and superfluid gap that we 
use are from T = QMC results J^iilii^ 

Referring to Fig.[3l if we consider chemical equilibrium, 
where 2^sf = M j^'' + Mj^^ i then at finite T we find that 
we have both a flux of t atoms entering the superfluid 
region and a flux of [ atoms leaving the superfluid region 
(or a flux of J, holes flowing in the opposite direction). 
However, atoms, not holes, are removed during evapora- 
tive cooling and thus we require the flux — across the 

interface. This is achieved by setting 2^sf ^ /^j-^^ +/^|^^ 
such that the concomitant mass current flowing into the 
superfluid cancels the J, hole current, so there is only a 
net motion of t atoms. Thus, the superfluid core ensures 
that the flux of evaporating atoms remains spin-polarized 
even when (/Lt| — /J-x ) is suppressed at the trap center due 
to local evaporative depolarization. 
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V. CONCLUDING REMARKS 

To summarize, we have shown that evaporative cool- 
ing of a trapped, partially-polarized, high-aspect-ratio 
atomic Fermi gas can create a spin current which in turn 
can produce the large 5c seen in the Rice experiment. 
Moreover, when this low-temperature gas cloud develops 
a superfluid core, we find that the atom flux across the 
superfluid-normal interface remains spin-polarized, with 
a small drop in the average chemical potential across the 
interface. An alternative equilibrium explanation for the 
large 5c in the Rice experiment is that finite-size effects 
in a high-aspect-ratio trapped gas are responsiblej ^^'^^i^^ 
However, microscopic studies of the surface tension at 
the superfiuid/normal interface in a trapped gas suggest 
that this scenario is not quantitatively consistent with 
experiments)^ Two possible avenues for the experiments 
to differentiate between equilibrium and non-equilibrium 
scenarios are: (1) To explicitly look for relaxation (on 
time scales that we estimate to be < 0.1 sec) towards 
equilibrium after deepening the optical trap and thus 
slowing or stopping the evaporative depolarization; and 
(2) To move the trap lip off-center with respect to the 
optical trap minimum by moving the minimum of the 
magnetic potential. If our non-equilibrium explanation 
is correct, at the lowest temperatures the superfiuid core 
should form where the local evaporative depolarization 
is occurring at the trap lip; at equilibrium the super- 
fluid will instead form at the minimum of the overall 
trapping potential where the density is highest. Indeed, 
a very recent experiments^ by Salomon's group at the 
ENS already casts some doubt on the alternative equi- 
librium explanation: their trap aspect ratio ('^ 22 : 1) 
approaches that used in the Rice experiment, yet they 



observe a 5c that agrees with the MIT result. Whilst it 
is still possible that this lower 5c may be because the tem- 
perature in the ENS experiment is higher, another key 
difference between the Rice and ENS experiments is that 
the time scale for evaporation in the ENS experiment is 
longer,-- a fact which appears to be consistent with our 
non-equilibrium proposal. 

More generally, the non-equilibrium scenario we have 
analyzed here suggests a route towards the realiza- 
tion and investigation of spin currents in a strongly- 
interacting Fermi gas. Moreover, it illustrates the impor- 
tance of non-equilibrium phenomena in cold-atom sys- 
tems. In typical experiments with ultracold atoms, the 
dynamic range between the microscopic time scales for 
atom motion, atom-atom interactions and the total du- 
ration of the experiments is not so large, often 5 orders 
of magnitude or less. This means that full equilibration 
of a large cloud of atoms might not be possible under 
some circumstances. Here we have discussed a trapped 
unitary Fermi gas, where the microscopic collision rate at 
low temperature is roughly sp jh. When one instead con- 
siders atoms in an optical lattice, the microscopic rates 
for atom hopping or for superexchange can be substan- 
tially slower than this, making it even easier to drive the 
system out of equilibrium.-^^ 
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